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I. INTRODUCTION 



The Haldane-Shastry spin chain ||T],0 is a one- dimensional integrable lattice model with 
interaction which has provided a number of valuable notions in low-dimensional 
quantum systems. Among others, (fractional) exclusion statistics proposed by Haldane 
has received considerable attention. From the viewpoint of exclusion statistics, the Haldane- 
Shastry model is regarded as an idealized model which is completely free from any irrelevant 
perturbations and describes the fixed point Hamiltonian for massless spin chain models. In 
this idealized circumstance, it may be possible to observe characteristic features inherent in 
exclusion statistics. 

Motivated by this, we have recently extended the Haldane-Shastry model to include 
twisted boundary conditions, and have shown that exclusion statistics can be explicitly ob- 
served in the period of the spectral flow ||^. The strategy to probe exclusion statistics has 
been based on the idea that the response to twisted boundary conditions (or equivalently 
to external gauge fields) enables us to extract the knowledge on various properties for in- 
teracting particle systems ||7|-[TT|. As mentioned above, the Haldane-Shastry model is free 
from any irrelevant perturbations, allowing us to observe exclusion statistics clearly in the 
spectral flow. It has been indeed shown P that the period of the ground state is solely 
determined by the statistical interaction in exclusion statistics. 

In this paper, we generalize the twisted Haldane-Shastry model to the case with SU(z/) 
spin symmetry, and solve it exactly. Such a generalization to the multicomponent models 
may be interesting in its own right [0,0 and also its close relationship to the fractional 
quantum Hall effect. We show that the period of the spectral flow in the SU(z/) model is 
indeed determined by the matrix of the statistical interaction in exclusion statistics, i.e. the 
SU(i/) Cartan matrix. 

In the next section, we introduce the Hamiltonian of the SU(z/) Haldane-Shastry model 
with twisted boundary conditions. We explain how the twisted boundary condition is im- 
posed consistently with the long-range nature of the interaction. In section HI, the 
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exact solution of the twisted SU(i/) model is derived by exploiting the Jastrow-type ansatz. 
In section IV, the exact spectrum thus obtained is shown to be reproduced by the Bethe 
ansatz. We discuss the spectral flow of the energy spectrum as a function of the twist angle 



in section V. By resorting to the notion of the motif in the Bethe ansatz, we give an 
interpretation for the period of the spectral flow in terms of exclusion statistics. Summary 
and discussions are given in section VI. 

II. MODEL HAMILTONIAN 

The Haldane-Shastry model has long-range interaction, so that it is not trivial to 
impose twisted boundary conditions without loss of integrability. Recently, we have proposed 
one of the ideas for this purpose p. We flrst recall that the Haldane-Shastry model [jl],|^ 
was originally introduced by imposing the periodic boundary condition and by summing up 
the pair- wise interaction around the ring inflnite times From this point of view, we 
start with the following Hamiltonian, 

oo 

H = Z-^ 1 'i T\j\1 7i'Pn,n+mN 

rf^n'mt^oo [u - u' - mNf 2 

oo -| 

= E E 

n<n' m=-oo 



{n — n' — mNY 



X 



^ QG*+ \k=\ ^'^Z 



where n denotes the index for sites n = 1,2, ■■■,A^ and Pn-n' = '^HaTnT^' + l/'^ is the 
exchange operator in SU(z/) algebra. In the second line, is the set of the positive roots 
of SU(z/) algebra, the generator E'^ is the so-called step operator associated with the root a 
at site n, and with = 1, 2, ■ ■ ■ , i/ — 1 is the hermitian generator in Cartan subalgebra. 
We flx the normalization of the generators as tr(T'^T'') = \5ah- The Hamiltonian (2J.) was 



already solved in the case of the periodic boundary condition |12,13|, where the summation 
over m can be explicitly carried out to give the well-known sine-inverse-square interaction 
H. However, once we impose twisted boundary conditions, it becomes a non-trivial solvable 
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model compatible with twisted boundary conditions. The form of the Hamiltonian 
implies that the Cartan basis may be a natural basis to treat the system with twisted 
boundary conditions. By this reason, we shall span the basis of the wave function in terms 



of the Cartan basis, which is different from that in [^], though both bases give the same 
results in the case of the periodic boundary condition. 

To be more specific in our notations, let us denote the simple roots by with i = 
1, 2, ■ ■ ■ , 1/ — 1. Then the positive roots can be expressed as a^k) + ■ ■ ■ + with 1 < 
k < I < u — \. Step operators associated with the positive roots are denoted simply as 

= for A; = 1, ■ ■ ■ , z/ - 1 and E^^'^^ = ^^c^+d+ '+^o) for 1 < A; < / < z/ - 1. 

We label the state at each site as which is specified by the u weight vectors with 

i = 0, 1, ■ ■ ■ , 1/ — 1 in the fundamental representation. As we have fixed the normalization of 
the generators, the weight vectors satisfy 



I - ^ for z = J 



Hi) ■ m = { (2-2) 

for t ^ J. 

The highest weight state is assumed to be |/U(o)), and other descendant states are created by 
the lowering operators E^^ such that |yU(fc)) = -E'"^|/i(o))- Therefore, in this representation, 
E^ connects directly the highest weight state |/i(o)) with a descendant state \fi{k)), while 
]^(k,i) connects two descendant states \fi{k)) and |yU(«)). 

The boundary condition we impose on the Hamiltonian in this paper is 

rp±k _ ±2iTi(p^m rp±k 
^n+mN ~ ^ i 

Hn+mN = (2.3) 

where the first line is the definition of the v — 1 twist angles 0^ and the second line naturally 
follows from the commutation relation of the step operators. In what follows, we refer to the 
angle (pk defined in unit of 27r as the twist angle. Note that we can impose twisted boundary 
conditions independently on each species. 



In order to obtain the exact solution, it is convenient to introduce the periodic operators 
by the following gauge transformations 

ipifc ±2-ni4>^n/N rp±k 

In the remainder of the paper, we always use the step operators in the gauge transformed 
form. Now we restrict ourselves to rational values 

0fc = - (2.5) 

qt 

with pk and being integers. This constraint is, at present, essential to solve the twisted 
model. Then the summation over n can be done explicitly, and the Hamiltonian (|2.1|) is 
now written as 

H = {^)' \iT + H,^,l (2.6) 

where T and iJint are the hopping and interaction terms, respectively. The former is defined 
by 

T = Y.i:J^Sn- n')E'^E~f^ + J^.-^.i^ - n')i?f ^i^^/'^^ (2.7) 

k n^n' k<l nj^n' 



Here Jrf){n) is the coupling compatible with the boundary condition ( p.3| 



j^(n) = 1 e2-^^'("+'"^)/''^sin- 



1^ m=0 



TT{n + mN) 



qN 

For later convenience, we rewrite the above expression as follows: 



for = ^. (2.8) 



T = ^rfc(0,) + i^TfcK0H), (2.9) 

k ^ ki^l 



where Ski = (j)k - (J)i and 



^fc(0) = E - ^')EnE-'\ Tki{4>) = UiTkUl (2.10) 

Here Ui is the unitary operator of the vr rotation around + ■■■ + a(k)i i-e., Uk = 
e'x^{i'^{E^ + E~^)}. The interaction term is 



i^int = E Un - n') E H^nH'w + ■ (2.11) 

n^n' \k=l '^'^ / 

Note that Jo(^) is of the usual sine- inverse-square form Join) = sin^^ This completes the 
introduction of the twisted SU(z/) Haldane-Shastry model. We can see that after the gauge 



transformation (|2.4|) the effect of twisted boundary conditions is incorporated in effective 
hopping, and hence we can solve the Hamiltonian with periodic boundary conditions. 



III. EXACT SOLUTION 

In this section we construct the exact eigenstates for the model ( p.6| ) and obtain the 
corresponding eigenenergies. In the previous paper [Q, we have solved the simplest case, 
i.e. the twisted SU(2) model, by using the Jastrow-ansatz wave function. In that case, 
we have demonstrated that the Jastrow wave function can still be an eigenfunction of the 
twisted Haldane-Shastry model if we convert the operators into the periodic ones via the 
gauge transformation. In what follows, we will show that this is also the case for the SU(z/) 
model, and then obtain its exact solution. 

Before starting the calculation, let us fix some notations. There are u different states at 
each site. We take the highest weight state |yU(o)) as the reference state (background), and 
regard other descendant states \fi(k)) (fc = 1, 2, ■ ■ ■ z/ — 1) as particle states. Since the number 
of states in |/i(fc)) is conserved, let us denote it as M^. Then the relation J2'^zlMk = N 
holds. Let the subset of the sites denoted by {n^a^} with = 1, ■ ■ ■ , Mk be the positions of 
the state |/i(fc)). 

We now propose the following Jastrow-ansatz state as a candidate for the exact eigenstate 
of the twisted SU(z/) model, 

i^) = E E ^(■■■,K^},---;---,4,---)n n ^:(')i/^o,/^o,---,/^o), (s.i) 

fc=l„W<r...^r,''=' k=lak=l 
"l ^ Mk 

where ip is a. Jastrow-type wave function 

u-l Mk . u-l 

^ = n n n n - <v n n ^(-s - -S)- (3-2) 

k=lak=l k=lak<l3k k<loik,ai 

6 



Here, d{n) = sm{7rn/N), z = exp{27ri/N), and the current for each species is defined by 

N -Mk-Mo , , 

Jk = mod 1. (3.3) 

We wish to prove below that this is indeed an eigenstate of the Hamiltonian. 

Let us start by evaluating the action of the hopping terms on the wavefunction. We 
have introduced two kinds of the hopping Hamiltonians Tk{4>k) and Tki{(f)ki) in eg. ( p.lOp - 



The former exchanges the pairs of {n^^^} and {n^^^}, while the latter exchanges the pairs of 
{n^''^} and {n^'^}. In what follows, we first evaluate the action of Tk on the wave function, 
which can be calculated directly. Although the direct manipulation of Tki turns out to be 
difficult, there is a remarkable trick to simplify the calculation, found by Wang, Liu and 



Coleman [|15[ for the supersymmetric t-J model. By generalizing this trick to the SU(i/) 
model, we then evaluate the action of Tki by the use of the results obtained for the action 
of Tk. 

Although the wavefunction (|3.2| ) contains the coordinates of various species, we note that 
Tk acts only on |yU(o)) and |/U(fc))- We thus classify the set of total sites into three subsets 
{n^^^, {n^^^} and the remainder {n} — {n*^'^-'} — {n^'^^}. To simplify the expressions, we 
introduce the following notations for the elements of each subset, 

x^E{n^'^}, « = l,2,---,Mfc 

2/i e M - - {nW}, j = l,2,---,N-Mk-Mo. (3.4) 

Then, the problem reduces essentially to that of two species with the background, and one 
finds that the action of Tk on (|3.2|) is analogous to that of the spin-hopping term for the 



supersymmetric t-J model with interaction |jT6|. Therefore, by exploiting the techniques 
developed there, we have 

'^'^'^^'^ = E^0.H^'^"E n < n (3-5) 



where 



n=l a=l/3(^Q)=l 3=1 



_ 1 An) in) _ (1 - Z'^)zl + Z '')z]f 

B^p-l-g^„ 9^,- 1^^^2 (3-6) 



with Zn = z^" and 



= cos — + sin — cot e„, = ^^^^. (3.7) 



The last factor F in eq.( p.5| ) describes the interaction between different species. Note that 
such interaction terms automatically disappear in the case of SU(2). Expanding the above 

formula in power series of z"-' and 2;"", and using the original notation for sites, we finally 
find the expression for the action of T^, 

Il!M!t = j2wt\ (3.8) 

i=i 

where 

W^'^ = 2M,e{Jk + <Pk) + \Mk{N^ - 1) + \mu{MI - 1) + ^M^M^, (3.9) 

3 o z 

W^f = \ E ^o«) - <)) - \ E Un^'l - n'S) + ^ E E Mnl'! - n^), (3.10) 

= 2t ih + 0fc - y) E cot e„,.„„, (3.11) 

^ E cot e„^,/3^ cot e^^^ao + ^ cot 600/30 cot e«o„^. (3.12) 

The function £( J + 0) is 

eiJ + d)), for integer J, 

e(J + 0) =<! \ ^ (3.13) 

\ (e{J + - i) + £( J + + i) - i j , for half-integer J, 
where 5 (A;) is defined for every rational k as 

s{k) = {2[k] -N+l)k-[k] {[k] + 1) . (3.14) 

The expression ( p.8| ) holds for the currents which satisfy 

iV + Afc - M„ _ ^ ^ _^ ^ ^ _ / + A-4 - „ _ X ^ (3.15) 



2 - " V 2 

The derivation of the above equations is outlined in Appendix. 

We now move to the action of T^. As already mentioned, Tki exchanges the pairs of {n^'^)} 
and {n*-'^}, making its action on the wave function ( |3.2| ) rather complicated. However, we 
note the following identity 
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= A,Ij{- ■ ■ , {n(o)}, . . . ; - J, + iV/2, ■ ■ ■ , iV - 4, ■ ■ ■ , J,_i - Jfc + N/2), (3.16) 

where A is a constant independent of coordinates. Let us denote the r.h.s. of the above 
equation as Aipij^y Combining these two relations, we can evaluate the action of T^i quite 
similarly to that of T^. The key formula is 

which was previously found for the supersymmetric t-J model [|l^. The relation implies 
that what we have to do is to substitute the following replacements to eq. 



,n I Jkl + T foik^l 

{^(0)}^|^W} (Mo^MO, Jk^\ ' (3.18) 

^N-Jk ioTk = l, 

where Jki = Jk — Ji- The results are 

Td^ = j2wt\ (3.19) 

where 

W^f ) = 2Mke{Jki + N/2 + 0^,) + Im,{N^ - 1) + ^Mfc(M2 - 1) + ^M^Mf , (3.20) 

3d/ 

W^f ^ = I E ^o(nW - n^S) - I E Jo«) - <)) + i E E ^o«) - n^), (3.21) 



^ = 2i (Jh + 0h) E cot e,,,,, (3.22) 



r{kl) 

— \^kl -r Vkl) 

p^^H) ^ E cot 6a^/3^ cot Ga^a, + E COt 9a;/3; COt Ga^a^ , (3.23) 



provided that the currents satisfy the condition 



Mk- Ml Mk- Ml , , 

^ - 1 < Jfci + < 1 -. (3.24) 



Consequently, from eqs.(3.S) and (3.19), we have the action of the total hopping Hamiltonian 



^ on the wave function. 
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^-j:W., (3.25) 

^ k=l 

where 

Wi = E W,^'^ + for t = l,2, 3, 4. (3.26) 

k=l k,l=l 

In this expression, there still remain unwanted two- and three-body terms. In what follows, 
we show how these unwanted terms indeed vanish. 

First, we can easily confirm that the two-body terms depending on the currents, 
vanish. The other two-body terms W2 are reduced to 

- E Un'^l - 4?) - Te E ^o«) - <)) 

k=l aky^l3k k=l ak,ao 



+1(--2)EE E ^o«)--) 



I E E -^o^) - 4^j) - ^MoiN' - 1) + ^(i. - 2) X: M,(iV2 - 1). (3.27) 



Next, the three-body terms W4 are calculated as 



W4 = ^ cot Ga,^, cot Oa^ai 
fcj=0 



= J2 cot ©afcA E cot Qa^n " ^ot 0„,^, 

= - E E Mn'£ - nf^) + i E M,(M| - 1). (3.28) 



Combining these formulae, we end up with 

^ = 2£;-i|;jo«l -»£'). (3.29) 

Now, it is seen that the two-body terms are exactly canceled out with i^int. Consequently, we 
can prove that is an eigenf unction of the Hamiltonian, and the corresponding eigenenergy, 
defined by Hil) = (^-^j i?^, is given by 
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u—1 I/—1 

E=Y. MkiiJk + 0fc) + ^ E MkiiJki + N/2 + (Pki) 

k=l ^ k,l=l 

+ ^('^ + 4) E Mk{Ml + l -Mk- M,)Ml 

k=l k=l 

+ ^(3z/ - 2)(iV - Mo)(iV2 - 1) + l(iV - Mo)M2 + \mo{MI - 1) - 4^o(iV' - 1), (3.30) 

where the energy depends exphcitly on {0} via the first hne of the expression. This is the 
exact eigenenergy for the SU(i/) twisted Haldane-Shastry model, which is the main result in 
this section. 

It is now instructive to write down the ground state energy. If we restrict ourselves to 
= vM with an even integer M and < 01 < 02 < ■ • ■ < ipv-i < 1, the ground state is 

realized by the choice Mq = Mi = ■ ■ ■ = M^^i and Ji = J2 = ■ ■ ■ = Ju-i- Then the ground 

state energy as a function of the twisted angle reads 

^g.s. = M E - z/ + l)0fc - ^ly^iy -2)M^ - ^v{2v - l)M. (3.31) 

k=i 

Note that this is the ground state energy for restricted twist angles, from which it is difficult 
to get full information on the the spectral flow beyond the restriction, though it is possible 
in principle. The main difficulty comes from the expression of the spectrum ( p.30| ) which 
contains the Gauss symbol (see the definition of eilz) in ( ^.14| )). Namely, whenever one of 
{0fc} jumps over an integer or other 0j's, the energy has a singular cusp structure (see also 
Fig.l). It is hence not easy to pass through these cusps and trace the spectral flow. In order 
to overcome the difficulty and to discuss the spectral flow correctly in the full region of the 
twist angles, it is necessary to resort to an alternative method, i.e. the Bethe ansatz, which 
will be described in the following section. 

IV. BETHE ANSATZ APPROACH 

Before discussing the spectral flow, we flrst show that the exact spectrum obtained in the 
previous section can be reproduced by the Bethe ansatz (sometimes called the asymptotic 
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Bethe ansatz [Q). The advantage to exploit the Bethe-ansatz description is, as will be shown 
below, that the spectral flow is naturally determined by appropriately choosing the quantum 
numbers in the Bethe equations. Also, by using the motif picture which can visualize 
how the quantum numbers in the Bethe ansatz are determined, we can naturally interpret 
the period of the spectral flow in terms of fractional exclusion statistics. 

In order to formulate the Bethe equations for the SU(z/) model, let us introduce the u — l 
kinds of rapidities, denoted here as k^^^ with z = 1, 2, ■ ■ ■ , i/ — 1, and = 1, 2, ■ ■ ■ , J2a=l ^a- 



Following a standard procedure in the nested Bethe ansatz , we can formally write down 
the Bethe equations as 

^ "2 bi{^ai) 



\ E sgn(A;2)-A;«) + /l? + 0,-0._, 



2 



\ E ^Mk^ - k^Z^) + ^ E sgn(A;« - A;2:^ 



I E sgn(fc(-^i) - kt_:^) + Itl^ + - = I E ^Mktl^ - t?) (4.1) 

with i = 2, ■ ■ ■ , z/ — 2. Here we note that the two-body phase shift takes the step-like form, 
which is inherent in the systems |Q . This special form of the phase shift is essential for 
the model to be ideal in the sense of exclusion statistics 0. The total energy of the model, 
defined in the unit of (^-fj^ as before, is given by 

. A/i+--+M,^_i 

E = -N{N'-1)+ Y: ^(^i'.^) (4-2) 

ai=l 

with e{k) defined in ( p.l4| ). The energy thus seems to be solely determined by the rapidity 
k^^\ but we should remember that its configuration is affected by the other rapidities via 
the nested equations. Therefore, when we follow the spectrum as a function of the twist 
angle, we must specify the precise configuration including k^'^\ ■ ■ ■ , k^'^~^\ whose behavior 
should be quite different from each other. 
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To see clearly how the rapidities are nested with each other, we take the simplest but non- 
trivial example, i.e. the SU(3) case. It is straightforward, though a little bit complicated, 
to generalize the following discussions to the SU(z/) models. To be specific, we use the 
simplified notations, k^^ = ka and k^^ = Aj, and write down the Bethe equations as 

ka = /^'^ + 01 - J E sgn(/ca - + ^ E sgn(A;a - A,), (4.3) 
I E sgn(A. - A,) + Jf + 021 = lj:^gn{X, - k^), (4.4) 

where 02i = 02 - 0i- 

Here we notice that eq.( [4.4| ), which has been formally deduced by the nested Bethe 
ansatz, does not seem to hold for a fractional value of 02i at a first glance. This problem 
stems from the fact that the two-body phase shift in the present system is of the step-like 
form. We now wish to explain how we can deal with eq. (|4.4|) correctly. To this end, let us 
consider the following configuration at 0i = 02 = as an example, 

• • ■ < ^Q-i < Aj < A;„ < • • • . (4.5) 

If we put 021 = 1, we see from eq. (|4.4|) that if^ is shifted to if^ -|- 1 and the above 
configuration should be changed to 

• • ■ < /cq,_i <ka<\i<--- . (4.6) 

Namely, the position of Aj is exchanged with that of k^ sitting on its right neighbor. The 
fractional 02i between and 1 should smoothly interpolate these two configurations. One 
readily notices that this is done by introducing an infinitesimal width ri in the step-like 
phase shift and then taking the limit of ^ 0. This observation naturally leads us to divide 
the l.h.s. of eq. (|4.4|) as if^ -|- 02i = {if^ -|- [02i]) + (02i — [02i])- Then the fractional part 
021 — [02i] can be absorbed into sgn(Aj — k^) = sgn(O) in the case of Aj = ka for < 02i < 1. 
Consequently, we get the simple result from eq.( [4.4| ) 



^sgn(Ai - ka) = 021 - ^ (4.7) 
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for < 021 < 1. By substituting the above formula into ( [4. 3D , we finally have 

+ for ■ ■ ■ < ka-l < ka < ■ ■ ■ , , ^ 

ka = i (4.8) 
[koa + (f)2, for ■ ■ ■ < Xi < ka < ■ ■ ■ , 

for < 021 < 1, where fcoo is the rapidity for {</)} = 0. From this equation we can see that k 

with (without) A in its left neighbor labels the rapidities for species-1 (-2) degree of freedom. 

Up to now, we have restricted ourselves to the case 02i ^ 0. For the other case 02i < 0, 

A exchanges the position with k sitting on its left when 02i decreases. Namely, the role of 

"left" and "right" in the previous discussion exchanges. 

The above interpretation of the two-body phase shift can be generalized to the SU(i/) 

case. As a result, it is easily confirmed that by suitably taking the quantum numbers the 

Bethe equations ( [4.1| ) indeed reproduce the exact energy ( |3.3CI| ) for the SU(i/) case obtained 

in the previous section. 

V. SPECTRAL FLOW AND EXCLUSION STATISTICS 

We now discuss the spectral flow with the aid of the Bethe ansatz description. For this 
purpose, the notion of the motif is particularly useful, allowing us to trace the spectral flow 
correctly. Let us first recall that the solution of the Bethe equations ( |4.1| ) can be specified 
graphically by the sequence of and 1 denoting respectively empty and occupied states of the 
momentum kj 0. Note that the unoccupied momentum, 0, is introduced to represent the 
repulsion effect of the two-body phase shift. This is the essence of the motif. For example, 
as was shown in ref. the motif for the SU(2) case reads 0101 ■ ■ -010, which implies that 
the two-body phase shift enlarges the spacing of rapidities twice as large as the free fermion 
case. For the SU(A^) case, one can find from the Bethe equations (4J) that the rapidity 



kj^^ is subject to the constraint that any configurations with more than z/ — 1 consecutive 
I's are prohibited. For example, the ground state of the SU(3) chain described by eqs.(4.3) 
and (|4.4| ) is characterized by the motif for ka 

011011 ■■■0110, (5.1) 
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where 1 denotes the occupied ka- In order to include the solution of the auxiliary rapidity 
Aj, the motif should read more precisely 

01ol01ol---01olO, (5.2) 

where we have introduced o which denotes the position of the occupied A,. We need this 
type of the motif when we discuss the spectral flow for the SU(i/) case {u > 3) with several 
different twist angles (pi. 

To see how well the motif picture works for our problem, let us discuss the spectral 
flow by taking the SU(3) model as an example. In Fig.l, we have shown the numerical 
diagonalization results of the spectrum for the finite system. These numerical results may 
be complementary to the exact results obtained in section III, since the latter provides a 
specific series of the eigenstates. We can see several characteristic properties in the spectral 
flow, e.g. the linear 0-dependence on the twist angel, cusp structures, etc. These behaviors 
seem rather peculiar, because the smooth 0^ dependence for small (f) is expected for ordinary 
spin models with short-range interaction. We find that the above characteristic behaviors 
are closely related to high symmetry of the system which is essential for the model to 
be an idealized model without any irrelevant perturbations. We shall discuss later that the 
linear 0-dependence indeed results from high symmetry of the model. 

Now, let us analyze the spectral flow by exploiting the motif picture. As mentioned 
above the motif for the SU(3) singlet ground state is given by 0110110 • • • 0110. When the 
boundary is twisted as 0i = 02 = 4>, the motif develops with the increase of the twist angle 
as 

0110110 • • • 0110 ^ 1011011 • • • 1010 ^ 1101101 • • • 1100 ^ 0110110 • • • 0110, (5.3) 

where we have omitted o for simplicity, because I's move in parallel, and o is always sand- 
wiched by the same I's. Here, the arrow means that unit is added; Scf) — 1. The corre- 
sponding spectral flow in Fig.l is: (a) — >■ (b) — > (c). From this interpretation with the motif, 
we can say that the period of the above spectral flow is 3. There may be another type of the 
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spectral flow, since there are two independent twist angles 0i and 02- For example, when 
01 = and 02 = 05 the motif changes like, 

OlolOlolO ■ ■ ■ OlolO OlOlolOl ■ ■ ■ lolOlo = lolOlolOl ■ ■ • lolOO 
^ lOlolOlol ■ ■ • lOloO = olOlolOl ■ • -olOlO OlolOlolO • ■ ■ OlolO, (5.4) 

as the twist angle is increased. The corresponding spectral flow thus has the period 3. In 
both cases, we can see that the period of the spectral flow in the SU(3) model is 3. Similarly 
we can analyze the spectral flow for more general SU(A^) cases based on the motif. For 
example, the period of the ground state turns out to be u for 0i = 02 = ■ ■ ■ = (pu^i = 0- In 
this way, we can determine the period of the spectral flow in terms of the motif picture. 

We now wish to discuss how the above description is related to the notion of exclusion 
statistics. According to Haldane the statistical interaction g in exclusion statistics is 
deflned as 

where and A^^ are the number of hole- and particle-states in species /i, respectively. The 
cases g^i, = g5^^ with (7 = and g = 1 correspond, respectively, to free bosons and fermions. 
A remarkable point is that the statistical interaction g is uniquely determined by the two- 
body phase shift in the Bethe equations. For instance, in the one-component system such 
as the SU(2) model, the statistical interaction g is given by the two-body phase shift via 
the relation |((yf — 1) ^S^iki — kj). Therefore, we can deduce from the Bethe equations 

( [4.1|) for u = 2 that the statistical interaction for the SU(2) Haldane-Shastry model is g = 2. 
Since the effect of the two-body phase shift is represented schematically by the motif, one 
can see that the above analysis based on the motif is directly related to exclusion statistics. 
Namely, the period for the spectral flow is determined by the statistical interaction g. Q 



^If we include the Ising-anisotropy A = g(g — l)/2 for the SU(2) model, the period may become 



g, as was indeed shown for (7 = 4 in |10| 
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For the SU(z/) Haldane-Shastry model, it is known that the statistical interaction is given 
by the [u — 1) x [u — 1) Cartan matrix of SU(z/) algebra p!3| , 



9 



2-1 
-1 2-1 



V 



(5.6) 



7 



0-1 2-1 
0-1 2 

which directly reflects the structure of the nested Bethe equations ( [4 .11 ). We can easily 
generalize the above analysis to the SU(z/) model. For example, the period of the flow for 

is z/, which is related with 



ill 



(5.7) 



In this way, the period of the spectral flow in the SU(z/) model is determined by the matrix 
of the statistical interaction. Note that the matrix g corresponds to the spin sector of the 



topological-order matrix characterizing internal structure of the quantum Hall system [IS 



VI. SUMMARY AND DISCUSSIONS 

Summarizing, we have proposed the SU(z/) Haldane-Shastry model compatible with 
twisted boundary conditions, and have solved it exactly. The spectrum thus obtained can 
be correctly reproduced by the Bethe ansatz solution. We have then discussed the spectral 
flow of the model in terms of the motif in the Bethe ansatz, and found that the period of 
the ground state is z/, reflecting fractional exclusion statistics. 

Finally, we wish to make some comments on an unusual linear 0-dependence of the 
spectral flow. We show here that this behavior is related to high symmetry inherent in 
the present system. For this purpose let us flrst remember that the periodic Haldane- 
Shastry model has Yangian symmetry which is larger than SU(z/). Twisting the model breaks 
Yangian symmetry as well as SU(z/). However, a nontrivial conserved quantity associated 
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with Yangian still persists even for twisted systems of finite size , and we expect that this 
conserved quantity controls the linear dependence of the spectral flow. This feature can be 
clearly seen when we consider the conformal limit of the model. To be more specific, let 
us consider the SU(2) Haldane-Shastry model. For the periodic Haldane-Shastry model, we 



can take conformal limit explicitly as follows |]T^ 



H^ns = Y.'^--J-nJn-: (6.1) 



ra>0 



where is the current operator of the k = 1 SU(2) WZW model. The level- 1 Yangian 



generators are also given by 



Ql = Y. e'^'^ : J'-J: ■■ ■ (6.2) 

n>0 



For the twisted model, a similar calculation for the Hamiltonian ( |2.6| ) leads to 

Htus = Hpns + <PQl (6.3) 

From this equation, we can clearly see that Q\ remains as a conserved quantity, reflecting 
high symmetry of the system, even when the boundary is twisted 0. We can also see from 
( |6.3| ) that the existence of this conserved quantity explains why we have encountered the 
linear 0-dependence in the spectral flow. 

We have been concerned with the model of rational twist angles in this paper. Within 
the present framework, we may apply the continued fraction approximation to discuss the 
case of irrational twist angles. It remains still an open question to construct the exact 
solution directly for the irrational cases, which is now under consideration. Also, the proof 
for integrability of the present model remains as an interesting issue to be solved in the 
future work. 
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APPENDIX A: 



In this appendix, we derive ( p. 81 ) in a bit detail. We concentrate on the case where 
Qk = N — Mk — Mo is an even number and hence the current Jk defined in eq.( p.3|) is an 
integer. It is straightforward to apply a similar calculation to the other case {Qk = odd 
number). 

The key formula is 



^ N-l 
^ n=l 



2 ' 



for s + t = 2 

for s + t = 1 

for s = t = 

for others, 



(Al) 



2£(J + 0) + i(Ar2 


provided t<J + (f)<N — s, which was first introduced in |10|, by extending the original 
one in Now let us expand eq.( |3.(j| ) in power series of 1 — and 1 — It consists 
of finite polynomials of the form (1 — 2;")*(1 — 2;^")* with < s + t < + — 1. 

Therefore, provided that Qk/2 + Mk ~l<Jk + (pk^N — ((5fc/2 + M^ — 1), terms with 
3 < s + t should vanish in eg. ( p. 5]) . We can also find that such terms in the second order as 



:i-z 



n\2 



(1 



or {(1 — z"') + {1 — z ")}^ vanish, so that it is sufficient to consider 



in the expansion 

E n 

a /3(^a)=l 



N-Mk-Mo 

n 



Mk 



aj 



/ . ya/3 + 9 i/a/3 i/a7 



^)} 5: cote 



+ iE(i- E 9^:lm 



z^) + (1 - 2"'')} E cot e„j- cot e„jv 



(A2) 



Applying the formula (|A1|) , we then have 
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^^^^ = 2M,£( Jfc + + \Mk{N^ - 1) + ^MkiMl - 1) + Ig.Mfc 

o o 2 

-2t{Jk + (j)k-N/2)J2cote^j 

— 2 ^ cot Qaf3 cot Qaj ~ 2 51 "^^^ '^^^ (^"^) 



So far we have used the simphfied notations defined in eq.(O). The remaining task 
is to convert the expression into the original notation. For example, in the third line, 
Y.a,j cot Qakj = Y^ok Sr=/ cot Qa^ar Consequently we end up with the results shown in 
the text. 
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Note added in proofs 

After submitted the paper, we have learned that Liu and Wang ( |cond-mat / 96080261) 
solved the supersymmetric t-J model with 1/r^ interaction for irrational twist angles. Their 
method can be directly applied to the present model. 
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FIGURES 




FIG. 1. Exact spectral flow of the iV = 6 SU(3) model with Mq = Mi = M2 = 2 and 
(f) = 4>i = 4>2- Lower 20 levels are plotted. The flow of the ground state (0110110) is described by 
(a) ^ (6) - (c). 
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